PROCEEDINGS 

OF SCIENCE 



HU-EP-07/32 



o 
o 

(N 

< 



Fermionic fields in the pseudoparticle approach 



Oh! 



> 
On 

in 
m 
CN 

od 
O 

o 



Marc Wagner* 

Humboldt-Universitat zu Berlin, Institut fiir Physik, Newtonstrafte 15, D- 12489 Berlin, Germany 
E-mail: }ncwagner@physik . hu-berlin . de 



The pseudoparticle approach is a numerical method to compute path integrals without discretizing 
spacetime. The basic idea is to consider only those field configurations, which can be represented 
as a linear superposition of a small number of localized building blocks (pseudoparticles), and to 
replace the functional integration by an integration over the pseudoparticle degrees of freedom. In 
previous papers we have successfully applied the pseudoparticle approach to SU(2) Yang-Mills 
theory. In this work we discuss the inclusion of fermionic fields in the pseudoparticle approach. 
To test our method, we compute the phase diagram of the 1+1-dimensional Gross-Neveu model 
in the large-rV limit as well as the chiral condensate in the crystal phase. 
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1. Introduction 

Recently there have been several papers proposing models for SU(2) Yang-Mills theory with 
a small number of physically relevant degrees of freedom. These models include ensembles of 
regular gauge instantons and merons [jl|, |2|, the pseudoparticle approach fH [| §], superpositions 
of calorons with non-trivial holonomy [||, |7|] and an ensemble of dyons The common basic 
principle is to restrict the Yang-Mills path integral to those gauge field configurations, which can 
be represented as a linear superposition of a small number of localized building blocks (pseudopar- 
ticles), e.g. instantons, merons, akyrons, calorons or dyons. 

These models have been quite successful, when dealing with problems related to confinement. 
First of all, the potential of two static charges is essentially linear within phenomenologically rel- 
evant distances. Moreover, a confinement-deconfinement phase transition can be modeled, and 
numerical results for various quantities, e.g. the string tension, the topological susceptibility, the 
critical temperature or the low lying glueball spectrum, are in qualitative agreement with results 
from lattice calculations. 

However, all these models exclusively consider pure Yang-Mills theory. Therefore, incorporat- 
ing fermions is an interesting issue. In this paper we present first steps in this direction: we propose 
a method how to deal with fermionic fields in the pseudoparticle approach, and we test this method 
by applying it to a simple interacting fermionic theory, the 1+1 -dimensional Gross-Neveu model 
in the large-A^-limit. 

2. Fermionic fields in the pseudoparticle approach 
2.1 Basic principle 

The starting point is action and partition function of any theory with quadratic fermion inter- 
action: 



where (j) denotes any type and number of bosonic fields, e.g. the non-Abelian gauge field in QCD, 
and Q is the Dirac operator, which, of course, depends on these bosonic fields. 

To stay close to the spirit of the pseudoparticle approach, we consider fermionic field configu- 
rations y, which can be represented as a linear superposition of a fixed number of pseudoparticles: 



Each pseudoparticle is a product of a Grassmann valued spinor r\j and a function Gj, which is 
localized in space as well as in time (the term pseudoparticle refers to this localization). The 
integration over all fermionic field configurations is defined as the integration over the Grassmann 
valued spinors rjj: 




(2.1) 



(2.2) 




(2.3) 




(2.4) 
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Integrating out the fermions yields 

S cfte ,veW>] = |^ +1 xif(0)-ln(det((G ; |G|G / ))) (2.5) 
Z oc f Dtye-SmvAvM^ (2.6) 



where the "fermionic matrix" (Gf\Q\Gf) is the Dirac operator represented in the pseudoparticle 
basis. We will refer to this pseudoparticle regularization as Q-regularization, and we will shortly 
point out that this Q-regularization is not suited to produce physically meaningful results. 

In the case that det(<2) is real and positive, det(<2) = y / det(2 t 2). This suggests another 
pseudoparticle regularization: 

S cftelV e[0] = y n ^ +1 xif(0)-^ln(det((G J |G t G|G / })). (2.7) 
In the following section we will argue that this <2 ^-regularization has significant advantages over 



the <2-regularization ([2.5D. 

Note that using eigenfunctions of the Dirac operator as "pseudoparticles" yields the well 
known finite mode regularization [Toj ]. 

2.2 The <2-regularization versus the <2^<2-regularization 



The problem of the <2-regularization (2.5) is that applying the Dirac operator Q to one of 
the pseudoparticles Gy in general yields a function, which is partially outside the pseudoparticle 
function space span{G„}: 

QG f (x) = Y,a fk G k (x)+h f H f (x) (2.8) 

k 

with Hji normalized and Hf _L span{G„}. If \Y*k a j'kGk\ 3> \hf\, the situation is uncritical. How- 
ever, as soon as \Y,k a fkGk\ ~ \hj>\, serious problems arise: when computing the fermionic matrix 
elements (Gj\Q\Gf), a significant part of QGf is simply ignored, namely hfHf, because it is 
perpendicular to the pseudoparticle function space span{G„}. 

On the other hand, the <2^<2-regularization ( |2.7[ ) has the following advantage: both the left 
hand sides (Gj\Q? and the right hand sides Q\Gf) of the fermionic matrix elements (G j\Q^Q\Gf) 
might be (partially) outside to the pseudoparticle function space span{G„}, but they form the same 
function space span{<2G„}, in which their overlap is computed. Of course, the above problem of 
partially perpendicular left and right hand side function spaces does not exist anymore. 

For more elaborate arguments, especially why one can expect to obtain correct results from 
the 2^2-regularization, we refer to [|TT|]. 



3. Testing the method: the Gross-Neveu model in the pseudoparticle approach 
3.1 The 1+1-dimensional Gross-Neveu model in the large A^-limit 



As a testbed for our pseudoparticle method we use the Gross-Neveu model [12], which is 



a four fermion interacting theory with ,/V identical flavors. Action and partition function of the 
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1+1 -dimensional Gross-Neveu model are given by 



Z = J (^f[D\if {n) D\ff^e- s , 



(3-1) 
(3.2) 



where /J. is the chemical potential and g the dimensionless coupling constant. To get rid of the four 
fermion term, one usually introduces a scalar field a. Integrating out the fermions yields 



Z oc / Dae 



y ^effective 



(3.3) 
(3.4) 



with A = Ng 2 . 

In the following we consider the large-Af limit, in which the model can be solved analytically 
fli~3[ , p] , |15| |. This amounts to using an infinite number of flavors N, while A = Ng 2 is kept constant. 
Note that in the N — > °° limit only a single a-field configuration contributes to the partition function 



( p.4| ) minimizing the effective action. Note also that in the large-Af limit o is proportional to the 



chiral condensate, i.e. a 



-s 2 Lti? ( "V n) - 



3.2 Numerical results: the phase diagram and the chiral condensate 

From a technical point of view computations in the pseudoparticle approach are quite similar 
to those in lattice field theory. The number of pseudoparticles corresponds to the number of lattice 
sites, while the distance between neighboring pseudoparticles plays a role similar to the lattice 
spacing. The scale can be set by any dimensionful quantity and it can be changed by choosing a 
different value for the dimensionless coupling constant. For a recent lattice study of the Gross- 
Neveu model we refer to [Hql. 



For the following computations we apply the £2^<2-regularization (1.1). As pseudoparticles we 
use a large number of uniformly distributed hat functions, more precisely B -spline basis functions 



B-spline basis functions, degree 2 





Figure 1: B-spline basis functions in one and two dimensions. 
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of degree 2 (cf. e.g. [|17|]), which are shown in Figure [j]. There is one fermionic pseudoparticle 
per unit volume, the spatial extension of the periodic spacetime region is L\ = 144 and the tem- 
poral extension Lo varies, corresponding to different temperatures T = 1/Lq. The main reason for 
considering such pseudoparticles is that they yield a sensible set of field configurations: they form 
a piecewise polynomial basis of degree 2, i.e. any not too heavily oscillating field configuration 
can be approximated. Therefore, if the pseudoparticle method we have presented in Section § is 
a useful numerical technique, we can expect to reproduce correct Gross-Neveu results. In other 
words, B-spline basis functions are suitable pseudoparticles for testing our approach. 

At first we perform computations of the chiral condensate a at chemical potential /I = and 
temporal extension Lq = 8 for various values of the coupling constant A . As it is in lattice calcula- 
tions different values of A correspond to different physical extensions of the spacetime region and, 
therefore, to different values of the temperature. From these computations we determine that value 
of A, where a just vanishes: A cntlca i = 1.153. For all further computations we use A = A clitical . By 
doing this we have set the scale, since from now on Lo plays the role of inverse temperature such 
that Lo = 8 corresponds to the critical temperature of chiral symmetry breaking. 

After that, we perform a low temperature computation at Lq = 48 or equivalently T = 7^. itical /6, 
to obtain an approximation of the zero temperature value of the chiral condensate: Oq = 0.221. 




b) 



phase diagram for homogeneous o 
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Figure 2: a) g/gq as a function of /i/Gq and T /Gq (red dots) together with the analytically obtained phase 
boundary (blue line) and the tricritical point (fj, / Go,T / 'gq) = (0.608 ,0.318) separating first and second order 
phase transitions (black dot), b) Phase diagram for homogeneous chiral condensate (red dots: pseudoparticle 
results; green line: analytical result), c) Two sections trough the phase diagram showing g/gq as a function 
of pi /gq at T /gq = 0.283 (first order phase transition) and T /gq = 0.377 (second order phase transition). 
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This allows us to express all dimensionful quantities in terms of Gq. 

Now we are in a position to compute the chiral condensate at arbitrary temperature T / Go and 
chemical potential /i/do- Results for homogeneous chiral condensate are shown in Figure |||a to- 
gether with the analytically obtained phase boundary [13, 14] and the tricritical point separating 
first and second order phase transitions. Pseudoparticle and analytical results are in excellent agree- 
ment both for the phase boundary (cf . also Figure ^|b) and for the order of the phase transition (cf . 
also Figure ||c, where we have plotted g/Go as a function of /J./ Go for two different values of T/go, 
one in the first order region and the other in the second order region). 

For inhomogeneous chiral condensate a third so called crystal phase appears [15], where the 
minimum of the effective action (33) is not anymore given by a homogeneous chiral condensate 
G. In addition to the fermionic fields we also represent a in terms of B-spline pseudoparticles (for 
details cf. [11]). As before, the pseudoparticle phase diagram and the analytically obtained phase 
diagram are essentially indistinguishable (cf. Figure ||a). 

We have also compared the pseudoparticle chiral condensate and the analytically obtained 
chiral condensate at various points (}i / Go,T / Gq) inside the crystal phase; again, there is excellent 
agreement. Figure [3]b shows the emergence of a crystalline structure: the kink-antikink structure 
close to the left phase boundary changes to a sin-like behavior, when approaching the center of the 
crystal phase. 

Note that we have performed the same computations also with the naive Q-regularization. As 
expected the results are completely wrong, e.g. there is no chirally symmetric phase even in the 
simple case of homogeneous chiral condensate. One can easily show that this is inherent to the 
<2-regularization and not a problem of the number or the type of pseudoparticles applied [11]. 



a) 



phase diagram for inhomogeneous o - X = 1 .153, L-, = 144 





0.6 




0.5 






jo 


0.4 


h- 




<D 




3 

cc 


0.3 


CD 




Q. 




tem 


0.2 




0.1 










exact phase boundary 




pseudoparticle results » 




\ a = 










crystal 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

chemical potential n/a 



b) close to the left phase boundary 

a(x,)/a - X = 1.153, n/o = 0.636, T/c = 0.141, . 



1 

o 05 
5 
° -0.5 
-1 



inside the crystal phase 

a{x,)la - X= 1.153, |i/o = 0.640, T/o = 0.141, 



1 

0.5 

-0.5 
-1 




20 40 60 80 100 120 140 



o(x,)/o - \= 1.153, n/o = 0.638, T/o ■ 0.141, ... 




20 40 60 80 100 120 140 



o(x,)/o - X= 1.153, MAr = 0.636, T/c = 0.141, 



1 

0.5 



■0.5 
-1 



20 40 60 80 100 120 140 



20 40 60 80 100 120 140 



Figure 3: a) Phase diagram for inhomogeneous chiral condensate (red dots: pseudoparticle results; green 
line: analytical result), b) The pseudoparticle chiral condensate for T/ao = 0.141 and different values of 
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4. Summary and outlook 

We have proposed a method to incorporate fermionic fields in the pseudoparticle approach. 
While the naive <2-regularization is not suited to produce any useful results, the Q^-regularization 
has the potential to yield correct and physically meaningful results. The computation of the phase 
diagram of the Gross-Neveu model with the <2 + <2-regularization both for homogeneous and for 
inhomogeneous chiral condensate has been a first successful test of the pseudoparticle approach 
applied to fermionic theories. 

The next step is to apply the pseudoparticle approach to QCD and to identify a small number 
of physically relevant degrees of freedom, probably fermionic pseudoparticles, which are able to 
approximate typical low lying eigenmodes of the Dirac operator. The goal is to obtain a model 
with a small number of degrees of freedom, which exhibits both chiral symmetry breaking and a 
confinement deconfinement phase transition at the same time. 
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